Phonon dispersion and lifetimes in MgB2 
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We measure phonon dispersion and linewidth in a single crystal of MgB2 along the F-A, F-M 
and A-L directions using inelastic X-Ray scattering. We use Density Functional Fheory to compute 
the effect of both electron-phonon coupling and anharmonicity on the linewidth, obtaining excellent 
agreement with experiment. Anomalous broadening of the E2g phonon mode is found all along F-A. 
The dominant contribution to the linewidth is always the electron-phonon coupling. 
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The discovery of 39 K superconductivity in MgB2 ^ 
has led to in-depth study of the material and a picture has 
emerged of a phonon- mediated superconductor with mul- 
tiple gaps [3, 0, , moderate electron-phonon coupling 
(EPC) SllEB: and anharmonicity ^,4,^1]. However 
no measurements exist concerning either phonon disper- 
sion or the evolution of phonon lifetimes over the Bril- 
louin Zone (BZ), due to the absence of large single crys- 
tals. Neutron scattering on powder samples [sl Il0l| has 
been limited to the determination of phonon density of 
states. Raman spectroscopy, which is restricted to the 
BZ center, has shown that the optical mode with E2g 
symmetry, corresponding to in-plane distortions of the B 
hexagons, is strongly damped 11. 12. .1,1. .14]. 

Phonon damping can be caused by (i) EPC medi- 
ated phonon decay into electron-hole pairs |l5j |. or 
phonon-phonon interaction due to anharmonicity 
The linewidth (the inverse of the lifetime) of a given 
phonon is the sum of both contributions. Direct deter- 
mination of the contribution of each phonon mode to 
EPC from the measured linewidth is on ly p ossible 
if the anharmonic contribution is negligible [ITj and is 
seemingly questionable for MgB2 where many calcula- 
tions IJjSH suggest strong anharmonic effects. 

In this work we present the first measured phonon 
dispersion curves and linewidths (where possible) along 
three directions in the BZ, F-A, F-M and A-L. We cir- 
cumvented the problem of sample size by using high reso- 
lution inelastic scattering of a focused and intense X-Ray 
beam at the European Synchrotron Radiation Facility 
(beamline ID28), a technique successfully used in sin- 
gle crystalline samples and in particular for the measure- 
ment of high energy optical modes To understand 
the mechanisms governing the measured phenomena, we 
calculated phonon dispersion, the contributions of EPC 
and anharmonicity to the linewidth and the structure 
factors using Density Functional Theory (DFT). 

Small single crystals of MgB2, suitable for inelastic X- 



ray scattering experiments have recently become avail- 
able. The crystal used in our experiment was grown at a 
pressure of 30-35 kbar. A mixture of Mg and B was put 
into a BN container in a cubic anvil device. The temper- 
ature was increased during one hour up to 1700 — 1800°C, 
kept stable for 1-3 hours and decreased during 1-2 hours. 
As a result plate-like MgB2 crystals were formed of which 
we used a sample of about 400 x 470 x 40 ^m? , with a 
measured in-plane mosaicity of 0.007°. The beam inci- 
dent on the sample was obtained from a high-resolution 
Silicon backscattcring monochromator using the (8 8 8) 
reflection at an incident energy of 15.816 keV. The X-ray 
beam was focused onto the sample by a toroidal mirror 
into a spot of 270 x 90 /im'^ (horizontal x vertical), full 
width at half maximum (FWHM). Slits before the sam- 
ple further limited the vertical beam size to 30 /zm. The 
scattered photons were analyzed in energy by five spheri- 
cal silicon crystal analyzers operating at the same reflec- 
tion order and mounted in pseudo Rowland circle geom- 
etry. The total energy resolution was 6.1 meV FWHM, 
as determined by a fit to a Lorcntzian lineshape. The 
momentum transfer Q was selected by rotating the 7 m 
long analyzer arm around the sample position, in the hor- 
izontal plane, which also contained the linear x-ray po- 
larization vector of the incident beam. The momentum 
resolution was set to 0.04A in the horizontal direction 
and O.O7A in the vertical direction. The following mea- 
surements were performed at a temperature of 300K: i) 
Q=(l 2 ^), in almost transverse configuration along the 
F-A direction, i.e. with Q ■ q « 0, q=(0 ^) being the 
phonon wavevector, ii) Q=(l-(5 1+^ 0), while nearly fol- 
lowing the F-M direction (0 < 5 < 0.05). iii) Q=(l-(5i 
2+^ 0.5+(52), while nearly following the A-L direction 
(0 < 81,82 < 0.04). The choice of BZ points measured 
was dependent on a series of conditions including the 
need to optimize the structure factors, single (F-A) or 
multi-analyzer (F-M, A-L) measurement mode and spec- 
trometer and time limitations. 
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FIG. 1: (color online) Energy loss scan in almost transverse 
geometry measured at Q=(l 2 0.3) corresponding to 0.6 F- 
A. The data, normalized to the incident flux, are shown with 
the least-squares flt (dashed line) and the ab-initio spectrum 
with and without broadening due to experiment and electron 
phonon coupling (solid lines). The broad peak corresponds to 
the damped E^g mode and is shown in greater detail in the 
inset. The peak at zero is due to diffuse scattering. 



In Fig. n we sliow data taken at the 0.6 F-A point in 
the BZ. The acoustic mode as well as the lower energy op- 
tical mode (Eiu) are visible as resolution-limited peaks. 
Most importantly, a broad peak is observed at higher en- 
ergy loss, corresponding to the i?2g mode. We performed 
least square fits to sums of Lorentzian functions with 
FWHM corresponding to the experimental resolution for 
the resolution limited peaks and a free parameter for the 
strongly damped phonon. These yield the dispersion as 
well as the linewidth variation over the BZ. Despite sta- 
tistical limitations (3-6 counts per minute on this peak 
along F-A) and tails of the peaks from the stronger, low 
energy phonons, the peak energy as well as the linewidth 
can be estimated with reasonable confidence. 

Fig. El shows a similar energy loss scan at (0.97 2.29 
0.54) close to the 0.58A-L point. A strong acoustic mode 
is seen at 30 meV. The peak at 50 meV corresponds to 
the branch and the one at 65 meV to another acous- 
tic branch. Finally two resolved features are seen at 85 
and 97 meV. These are the two E2g modes which in this 
region of reciprocal space are well separated from other 
modes. Though, given the statistics, it would be haz- 
ardous to estimate a linewidth, the comparison between 
the experimental and ab-initio spectra suggests that the 
linewidths of both the E2g modes are resolution limited 
and so the damping is much less than that along F-A. 
As for the measurement nearly along F-M, the structure 
factor for the optical modes strengthens only near the 
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FIG. 2: (color online) Energy loss scan, measured at Q=(0.97 
2.29 0.54) corresponding closely to 0.58 A-L, with the least- 
squares fit and the ab-initio spectrum calculated at 0.58 A-L. 

zone boundary. At the point measured nearest to M 
((1.05 1.45 0), not shown) the E2g and Eiu modes are 
comparable in intensity but only separated by about 1 
meV according to our calculation and we do measure a 
single peak only somewhat broader (FWHM « lOmeV) 
than the experimental resolution, indicating reduced i?2g 
linewidth. The proximity of the Eiu mode however pre- 
vents a firm conclusion in this regard. We mention that 
the calculated structure factors and energies show excel- 
lent quantitative agreement with our measured data of 
which we have shown only two examples. 

Similar analysis was done for several points along the 
three directions in order to experimentally determine the 
phonon dispersion and the linewidths. The difference in 
calculated phonon energies between the measured points 
and corresponding points exactly along F-M and A-L 
{6,61,62 = 0) is less than half a meV in all cases. We 
can thus compare the experimental phonon dispersion 
with the theoretical calculation along the high symmetry 
lines, as shown in the bottom panel of fig. 13 (circles). 

The measured intrinsic linewidth of the E2g branch, 
shown in the top panel of fig. 13 is strongly anisotropic in 
the BZ. Along F-A it is particularly large (ranging from 
20 to 28 meV), while near L and probably near M it is 
below the experimental resolution. 

Electronic structure calculations |2nj | were performed 
using DFT in the generalized gradient approximation 
[2l| . We used norm conserving pseudo-potentials [2^ . 
For Mg we used non-linear core corrections and 
we treated the 2s, 2p levels as core states. The wave- 
functions were expanded in plane waves using a 35 Ry 
cutoff. The calculations were performed with the exper- 
imental crystal structure, namely a = 3.083A and c/a = 
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FIG. 3: (color online) Bottom: Experimental (circles) and 
theoretical phonon dispersion (solid line) in MgB2 along F- 
A and F-M and A-L. In the region near the M point, the 
probable detection of the E2g mode is indicated with a square 
symbol (see text). The crosses indicate a parasite signal of 
unknown origin. Top: Intrinsic linewidth of the E2g mode. 
The experimental linewidth (circles) is large along F-A and 
below the experimental resolution (dashed line) near L and 
M. The theoretical result (diamonds) for the electron-phonon 
coupling contribution to the linewidth is also shown. Along 
A-L and F-M where the i?2g mode is non-degenerate both 
theoretical values are shown, when different. E2g linewidth 
decreases progressively from A to L and experimentally the 
two branches are resolved for the point nearest to L. 



1.142. The harmonic phonon frequencies were computed 
in the linear response 24] . We used a 16 x 16 x 16 
Monkhorst-Pack grid for the electronic BZ integration 
and first order Hcrmite-Gaussian smearing ^2^ of 0.025 
Ry. The dynamical matrix at a given point of the BZ 
was obtained from a Fourier interpolation of the dynam- 
ical matrices computed on a 6 x 6 x 4 phonon mesh. The 
resulting phonon frequencies are shown in fig. Hjand are 
in good agreement with a recent calculation [ll|. The 
agreement with experiment is remarkable. 

The contribution to the FWHM linewidth 7qjy at mo- 
mentum q for the i> phonon mode due to the electron- 
phonon interaction can be written as [l5| : 



l5kn,k-|-qmP'^('^kn)'^(£k+qm) (1) 



where the sum is extended over the BZ, Nk is the 
number of fc-points in the sum, and e^n are the en- 
ergy bands measured with respect to the Fermi level 
at point k. The matrix element is k+qm ~ 
{hn\5V / 6uq,y\h + qm) / y/2uj^ , where Uqi, is the ampli- 
tude of the displacement of the phonon v of wavevector 
q, LOqi, is the phonon frequency and V is the Kohn-Sham 
potential. 



In the calculations we used Nk = 30"^ inequivalent k- 
points and, in Eqs. 1^ and 0, we substituted the 6 
functions with Gaussians. The electron phonon coupling 
Aqi, is obtained from the linewidth as: 



(2) 



7V(o) = 0.354 states/ (MgB2 eV spin) being the density 
of states at the Fermi level. 

The second contribution to the linewidth is given by 
the anharmonicity in the crystal potential. At lowest 
order for the mode v oia. zone center phonon the FWHM 
linewidth is El S El: 
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duQi^dUqi^du^qri 



(3) 



E being the total energy. 



the Bose occupa- 



tion for mode /i at wavevector q, /q^^^ = (riq^ -r 
TT^qn + 1)<5 {ujQi, — ujqfi — LOqjj) dcscribes the decay in 
the two phonons /i and r], and Iq^^j^ = 2(rtq^ — 
nq^)6 (wojy — LOqf^ + ujqrj) dcscribes the 77-phonon absorp- 
tion and the /i-phonon emission. 

We computed the anharmonic linewidth at the high- 
symmetry points r. A, M. For the calculation at A we 
consider a 1 x 1 x 2 supercell with 6 atoms, while for the M 
point we use a 2 x 2 x 1 cell with 12 atoms. The third order 
matrices were evaluated using linear response theory and 
the 2ri -I- 1 theorem for metals |2^. The anharmonic 
contribution was evaluated at OK and 300i^. 

At F the anharmonic linewidth is largest for the i?2g 
mode and equal to 0.16 meV at T = OK and 1.21 meV 
at T ~ 300K. Both the values are negligible if compared 
with the experimental Raman linewidth of roughly 40 
meV |l3j , suggesting that the main source of broadening 
is the electron-phonon interaction. 

The results of the calculation of the two contributions 
to the linewidth at A and M are shown in Table|l| At the 
A point the linewidth of the E2g mode due to electron- 
phonon scattering is very large while the anharmonic con- 
tribution is more than an order of magnitude smaller. 
For sizable values of the EPC and at low temperature 
the anharmonicity is negligible showing that a measure- 
ment of the linewidth at q for the mode is essentially 
equivalent to determining Aajy. This is unexpected since 
earlier theoretical work 3, 4, 8] estimates anharmonicity 
to be important for a different but related quantity, the 
frequency shift. 

Since the anharmonic contribution is negligible at F, 
A, M, and computationally demanding, we only evalu- 
ated the EPC linewidth for the other points along the 
three directions. Using the calculated phonon frequen- 
cies, displacements and linewidths we compute the struc- 
ture factors for one phonon processes using the X-ray 
form factors, obtaining good agreement with experiment 
as shown in Figs. lll2l In the top panel of Fig. Owe show 
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samples, be measured for all modes over the whole Bril- 
louin zone so as to extract the anisotropic Eliashberg 
coupling function. 
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TABLE I: Calculated linewidths (meV) due to anharmonicity 
at OK (r2), 300K (F^"") and electron-phonon interaction 7^ at 
M and A for all modes. Phonon frequencies increase from top 
to bottom. Ai/ is the electron-phonon coupling. E-^g modes in 
boldface. 



q 


0.2 F-A 


0.5 F-A 


0.6 F-A 


0.8 F-A 


1.0 F-A 


\ expt. 


2.5±1.1 


2.6±0.6 


2.3±0.5 


3.6±0.7 


3.6±0.8 


\ theo. 


3.32 


2.80 


2.77 


3.12 


2.83 



TABLE II: Experimental and theoretical Aq^ of each of the 
two degenerate E2g modes along F-A. 



the theoretical results for the E^g branch. The result is 
consistent with the experimental value wherever the Eig 
branch is clearly visible namely along F-A and near the 
L point. Finally, the anharmonicity being negligible, we 
used Eq. ||2Jl to extract Aq^ of the mode along F-A 
using the measured linewidths and frequencies together 
with the calculated electronic density of states. In Table 
nil the experimental values are compared with the the- 
oretical predictions. The anomalously large EPC along 
F-A is due to the nesting factor of the B bonding p^.y 
Fermi surfaces, which are concentric cylinders centered 
on F-A 0, 111- The E2g modes, which modify the B-B 
distances, are the only ones with a sizeable matrix ele- 
ment, gj^jj k+qmi between electrons on these surfaces. 

In conclusion, we have measured phonon dispersion 
and linewidths in a sub-mm sized MgB2 crystal with in- 
elastic X-ray scattering confirming the power and ver- 
satility of this technique. Both acoustic and optical 
modes are detected and we find that the i?2g mode is 
anomalously broadened along F-A but that this broad- 
ening is not generalized over the Brillouin Zone. Our 
Density Functional Theory calculations of the dispersion 
and linewidth are in excellent agreement with experi- 
ment. They show that the dominant contribution to the 
broadening for all modes is the electron-phonon coupling, 
the anharmonic contribution being much smaller. Thus 
phonon linewidth in MgB2 is a direct measure of electron- 
phonon coupling and could, with the availability of larger 
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